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Currently, three types of superconducting magnets are used in particle accelerators: cos20, CCT, and serpen- 
tine. However, all three coil configurations have complex spatial geometries, which make magnet manufacturing 
and strain-sensitive superconductor applications difficult. Compared with the three existing quadrupole coils, 
the racetrack quadrupole coil has a simple shape and manufacturing process, but there have been few theoretical 
studies. In this paper, the two-dimensional and three-dimensional analytical expressions for the magnetic field 
in coil-dominated racetrack superconducting quadrupole magnets are presented. The analytical expressions of 
the field harmonics and gradient are fully resolved and depend only on the geometric parameters of the coil and 
current density. Then, a genetic algorithm is applied to obtain a solution for the coil geometry parameters with 
field harmonics on the order of 1074. Finally, considering the practical engineering needs of the accelerator 
interaction region, electromagnetic design examples of racetrack quadrupole magnets with high gradients, large 
apertures, and small apertures are described, and the application prospects of racetrack quadrupole coils are 


analyzed. 
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I. INTRODUCTION 


With the development of superconducting magnet technol- 
ogy in particle accelerators, cos20 coil and serpentine coil 
structures have been successfully applied to the final focus 
system in the interaction region. The cos20 magnet struc- 
ture is commonly used in the accelerator interaction regions. 
After decades of development, the cos20 superconducting 
magnets have matured for theoretical guidance and engineer- 
ing applications [1, 2]. The cos20 coil is saddle-shaped in 
the longitudinal direction and fan-shaped in the cross section 
around the beam aperture, which simulates the distribution of 
the quadrupole coil current J = Ipcos(20). The processing 
and manufacturing of cos20 magnets requires winding sys- 
tems, curing systems, assembly systems, and other related 
equipment. The magnet consists of coils, collars, and an iron 
yoke. The superconductor, spacer, and wedge were solidified 
as a solid whole after processing. Collars are used to pro- 
vide pre-stress and fix the coil to form a self-supporting struc- 
ture [3]. A small-aperture cos20 superconducting quadrupole 
magnet, QDO, is the basic scheme for the final focus system 
in the CEPC interaction region [4, 5]. A large-aperture cos20 
quadrupole magnet, MQXF, is composed of Nb3Sn super- 
conductors in the interaction region of the HL-LHC upgrade 
project [6, 7]. Serpentine winding, an innovation developed at 
BNL for direct-winding superconducting magnets, allows the 
winding of a coil layer of arbitrary multipolarity in one con- 
tinuous winding process and considerably simplifies magnet 
design and production [8, 9]. A winding machine with multi- 
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axis motion was used to wind each layer of the coil separately, 
and the four poles of each layer were wound simultaneously. 
Ultrasonic technology was used to heat the semicured glue 
and fix the superconductor in time. A low-temperature glue 
was brushed on the coil to hold it in place. After the low- 
temperature glue was cured, it was wrapped with a glass yarn 
on the outside to apply pre-stress, which eliminated the need 
for a extra pre-stress structure and reduced the magnet size. A 
serpentine quadrupole magnet was used for the BEPC-II Up- 
grade [10, 11], which was also proposed for manufacturing 
compact final focus magnets for the ILC. 


With improvements in the machining accuracy of machine 
tools, CCT magnets have also become the focus of accelera- 
tor magnet research. The position coordinates of each turn of 
the superconductor in the CCT coil are determined by a pa- 
rameter equation that is strictly arranged on the surface of the 
cylinder [12, 13]. The current density of the CCT quadrupole 
coil section also satisfies the cos20 distribution and the multi- 
layer coil skeleton structure can realize a combined functional 
magnet. According to the superconductor distribution, the 
CCT coil needs to slot the surface of the cylindrical skeleton 
using CNC machining or 3D printing technology, and then 
place the superconductors directly into the slot and fix it [14]. 
CCT quadrupole magnets are supported by a comprehensive 
theoretical foundation and have been used in a preliminary 
study of superconducting quadrupole magnets in the Fcc-ee 
interaction region [15]. 


With a further increase in the beam energy in a large 
particle collider, higher-gradient superconducting quadrupole 
magnets are required in the interaction region. A higher field 
gradient indicates that the superconductors need to carry more 
current. When the peak magnetic field in the coil exceeds 10 
T, the mature low-temperature superconductor NbTi reaches 
its theoretical limit. The current-carrying capacity of the su- 
perconductor Nb3Sn with a higher critical magnetic field also 


e4 decreases sharply. However, with an increase in the magnetic 
65 field, the critical current density of the high-temperature su- 
es perconductor decreases slightly, which is a better choice for 
67 realizing high-gradient superconducting quadrupole magnets. 
High-temperature superconductors are more strain-sensitive 
than low-temperature superconductors and are unsuitable for 
cos20, CCT, and serpentine coils with complex spatial ge- 
ometries. Therefore, a racetrack coil with a simple spa- 
tial shape and gentle bending is suitable for use in high- 
temperature superconductors. The above is an advantage of 
using racetrack coils for high-field magnets with a peak mag- 
netic field greater than 10 T, for some low-field supercon- 
ducting quadrupole magnets with a smaller peak magnetic 
field, considering the radiation damage and heating, a high- 
temperature superconductor is also a necessary choice to pro- 
vide sufficient operating margins to ensure operating stability. 
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Iron-dominated and coil-dominated superconducting mag- 
nets are the two basic types. Owing to the iron pole satura- 
tion phenomenon, iron-dominated racetrack quadrupole mag- 
nets have a low magnetic field intensity, and the magnetic 
field quality is difficult to control, which cannot meet the 
requirements of large particle accelerators [16-20]. From 
this perspective, a coil-dominated racetrack quadrupole mag- 
net is a reliable method for achieving an ultrahigh mag- 
netic field gradient. Recently, studies have been conducted 
on coil-dominated racetrack quadrupole magnets. The US 
LHC Accelerator Research Program has completed the de- 
sign and processing of racetrack superconducting magnets 
SQO1 and SQO2 using the superconductor Nb3Sn and con- 
ducted quench tests [21, 22]. However, the research on coil- 
dominated racetrack quadrupole magnets is limited and re- 
quires further improvement. First, current research on race- 
track quadrupole magnets is based on finite element soft- 
ware, which lacks the guidance of theoretical formulas such 
as the Cos20 coil [23, 24]. Second, return coils are added 
to form racetrack coils, which increase the size of the mag- 
net cross-section. Current distribution patterns are complex, 
causing difficulties in the design and manufacturing of mag- 
nets. [25, 26]. 


The focus of this study was a magnetic field numerical cal- 
culation algorithm for a coil-dominated racetrack quadrupole 
coil. First, the field harmonics and gradient analytical ex- 
pressions of the two-dimensional racetrack quadrupole coil 
as a function of the geometric parameters and current den- 
sity were derived based on the line current theory, coordinate 
transformation, and numerical integration. Second, a genetic 
algorithm is applied to obtain a solution for the coil geom- 
etry parameters with field harmonics of the order of 1074. 
Subsequently, the three-dimensional analytical expressions of 
the high-order field harmonics in the racetrack quadrupole 
coil were obtained using the discrete summation algorithm. 
Combined with the finite element software, the accuracy of 
the numerical calculation algorithm was verified. Finally, 
combined with the practical engineering design requirements 
of accelerator magnets, electromagnetic designs of racetrack 
quadrupole magnets with high gradients, large apertures, and 
small apertures are described, and the application prospects 
of racetrack quadrupole coils are analyzed. 
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Il. NUMERICAL CALCULATION ALGORITHM FOR 
TWO-DIMENSIONAL RACETRACK QUADRUPOLE COIL 
MAGNETIC FIELD 


122 


According to the complex formalism, the line current [ 
in the position z = xo + iyo generates a magnetic field 


127 B (z) = By (z) + iB, (z) in the position z = x + iy, which 
123 reads [27] 
Ipo 
B = — 1 
12 (a Qn (z — zo) ) 


where up is the permeability of the vacuum. Multipole field 
components in the aperture generated by line current J at po- 
sition zo = £o + iyo = roe’ are 


Tp1p Rrep 
By, = aeo ret cos no 
2m ro 
ti ae (2) 
I re . 
Åg = eal] f sin nO 
27 rọ 


14 6 Rye is the reference radius of the good region, deter- 
135 mined by the clear region of the beam. Normal and skew 
136 quadrupole magnets are the two basic quadrupole types. A 
137 normal quadrupole magnet was obtained by rotating a skew 
138 Quadrupole magnet by 45°. The multipole field components 
139 in the aperture generated by the rotating line current I’ at 
140 position 24 = zoe~*t = roe(%- 7) are 


Ino Bref T 

Pan sn(@, — — 

A In ri cos n(o p 
, Ea 6) 

Ip T 

Al a FO ref inns | 

n= Oe oe sin n(A9 7D) 


For the high-order field harmonics of the magnetic field, n 
in Eq. (3) equals 4N-2, N=/, 2, 3.... The Eq. (3) can be 
written as 


=i 
A’ = I Lo Reef 
n EA n 
T To 


T 


sin ((2n - D7 


(4N — 2)60) (4) 


When N is odd, n= 2, 10, 18... 


Ipo Rye 
Al = Ho ref oos ((4N — 2) 0o) 
2m ro 
147 I Rer (5) 
= — 0 ref og (np) 
2m ro 
= Bn 
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When N is even, that is, n = 6, 14, 22... 


nm—1 
E a Ino Reef 
no OR ro 


Ip Ry. a (6) 
“Po ref cos (nbo) 
27 


os ((4N — 2) 0o) 
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According to Eq. (5), and Eq. (6), 
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Br 
a Be = bn: m= 10,18... 
151 al, i (7) 
2 —" = bn = 6, 14,22... 
Bə ? n ? ? 


Among the high-order field harmonics ag, a@19 and a14, a 
153 negative sign must be added to harmonics ag and ay, in the 
154 transformation from a skew quadrupole magnet to a normal 
155 Quadrupole magnet [28]. 


152 
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Fig. 1. Schematic of current distribution in 0-45° quadrant of normal 
racetrack quadrupole coil. The three areas marked with different 
colors and numbers in the figure indicate that the definitive integral 
in rectangular coordinates needs to be solved in three regions. 


There are two ideas for the numerical calculation of a 
two-dimensional magnetic field: complex representation [29] 
and real representation. The real representation is also the 
method used in this study. Fig. 1 shows the current distri- 
bution diagram of a normal racetrack quadrupole coil in the 
quadrants from 0° to 45°. The coil cross-sectional distribu- 
tion makes the integration operation extremely difficult in 
coordinate systems. Although normal quadrupole magnets 
are commonly used in accelerator-focusing systems, a skew 
165 quadrupole magnet was considered as our research object to 
166 facilitate the integral calculation process in the cartesian coor- 
167 dinate system. When the skew quadrupole magnet is rotated 
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by 45° to form a normal quadrupole magnet, the multipole 
field also changes from the skew field A, to the normal field 
Bn, and such a change only affects the positive or negative 
signs of the field harmonics. Eq. (7) deduces the conver- 
sion relationship between the high-order field harmonics of 
the normal and skew quadrupole coils: 


Y 


d 


Fig. 2. Schematic of current distribution in 0-45° quadrant of skew 
racetrack quadrupole coil. 


Fig. 2 shows the current distribution diagram of the skew 
racetrack quadrupole coil in the quadrant ranging from 0° to 
45°. A parameter definition for the current block is presented. 
To derive the total field produced by the current block, the 
multipole field component expressions in polar coordinates 
were transformed into cartesian coordinates. The transforma- 
tion formulas for the two coordinate systems are as follows: 


raea 


0 = arctan (2) 


From Eq. (8), the multipole field components generated by 
the differential unit of the rectangular block are 


; RET 
HoJ 
dAn(block-abcd) =a ref nZ 


x sinn (arctan (*)) dydx 


(9) 


where j is the current density of the coil block. The mul- 
tipole field components in the aperture generated by the rect- 
angular current block shown in Fig. 1 can be calculated as 
-1 
Hoj Rre ref 


An(block— abcd) -f i On ‘Qn (a2 + yË 


x sinn (arctan (2 ) ) dyae 


(10) 


189 According to Eq. (10), each order multipole field A, of 
190 the entire racetrack quadrupole coil is obtained. For example, 
191 analytical expressions for the quadrupole field component and 
192 field gradient are obtained. Other analytical expressions of 
193 multipole fields Ag and Ajo are listed in the Appendix. 


Ay — Rreg p (PHP) (P+) 

o T (a? + c?) (b? + d?) an 
ga Az Moi (a? + d?) (b? +c?) 
Rrep m (a2 +2) (b2 +d?) 


15 The relative field harmonics normalized to the main 
1 quadrupole field component are defined as an = An/A2. 
197 Then, the high-order field harmonics ag, aio and a14 can be 
198 derived. These analytical expressions contain four parame- 
199 ters: a, b, c, and d. 

20 In the analytical expressions, parameter a represents the 
201 X-axis coordinate value corresponding to the left side of the 
202 rectangular current block, which is directly related to the 
203 aperture of the racetrack superconducting quadrupole coil. 
204 Therefore, after reserving the space required by the coil sup- 
205 port structure, the parameter a is the first determined variable. 
2% According to the aperture of the quadrupole coil in the CEPC 
207 interaction region, the parameter a is set to be 0.02 m. The 
20s reference radius R,-¢ was set as 0.0098 m [30, 31]. 

29 Subsequently, the relationship between the high-order 
210 field harmonics and other parameters was further studied. 
21 Among the high-order field harmonics of superconducting 
2 quadrupole magnets, the systematic harmonics ag, aio and 
213 a14 have been the main research objects. Field harmonics ag 
214 as functions of the coil parameters are shown in Fig. 3. The 
as different layered structures in the z direction are caused by the 
216 different parameters b. Different parameters c and d result in 
217 Significant changes in the surface of the high-order field har- 
218 monics bg. The intersection of the high-order field harmonics 
219 Surface and z=0 plane indicates that bg is equal to zero. At 
22 this time, the parameters b, c, and d corresponding to each 
2231 point on the intersection line are the solutions. 

22 Similarly, as shown in Fig. 3, the value range of the high- 
2233 order field harmonics ao exhibits a small change, and most 
224 results are within the design requirements. The intersection 
2235 Of the high-order field harmonics surface and z=0 plane indi- 
226 Cates that b10 is equal to zero. As shown in Fig. 3, harmonics 
227 a14 are extremely small over the entire solution area. The in- 
22 tersection of the high-order field harmonics surface and z=0 
229 plane indicates that b14 is equal to zero. 

230 From an analysis of the high-order field harmonics of ag, 
231 a10, and a44, it can be observed that aj, is relatively small. 
232 Thus, ag and ajo are the key to this study. Under the condi- 
233 tion that parameters a, b, and c are determined, the relation- 
234 Ship between parameter d and the high-order field harmonics 
235 is Studied. ag and ajo as functions of the coil parameters are 
23 Shown in Fig. 4. All the points in the solution set regions 1 
237 and 2 meet the requirement that the high-order field harmon- 
238 ics be less than 5 x 1074, which can be used as the solution 
233 Of the model. However, such a graphical solution is com- 
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Fig. 3. High-order field harmonics distribution of racetrack 
quadrupole coil 
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Fig. 4. Distribution of high-order field harmonics under the same 
parameters a, b, and c. 


240 plicated, and a coil pole composed of a single-block region 241 causes a large excitation current in the superconducting cable. 
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Therefore, the model must be further optimized to improve 
the efficiency of the solution and reduce the difficulty in mag- 
net design and manufacturing. Increase the current blocks in 
each pole and improve the ability to adjust the magnetic field 
quality. Meanwhile, the solution method was optimized, the 
efficiency of the model solution was enhanced, and the prac- 
tical application value of the racetrack coil analysis algorithm 
was improved. 

To make the design of the theoretical model more suitable 
for actual situations and increase the adjustment ability of the 
quadrupole coil, a double-layer quadrupole coil composed of 
two blocks was used in the theoretical model. A double-layer 
racetrack quadrupole coil is shown in Fig. 5. For the range 
0-45° in the first quadrant, the two rectangular blocks have 
the following constraints: 0 < cı < dı < a, < bı and 
0 < cg < dg < ag < bg. To ensure that the two rectangular 
blocks do not overlap, we add the following constraint ag > 
by: 


block a,b,c,d, 


block a,b,c,d, 


one pole 


Fig. 5. Racetrack quadrupole coil with two blocks in each pole. 


According to Eq. (10), the high-order field harmonics of 
the racetrack quadrupole coil consisting of two blocks were 
obtained: 


ü An —blocks 
n—blocks = GY 
A2- blocks 


(12) 
— An—block aıbıcıdı T An blockagbecods 


A2~block arbıcıdı + Å2—blockazbzc2d2 


According to the design requirements, the following ex- 
pression can be obtained: 


_ Ag—block airbıcıdı + Á6—blockazbzc2d2 


ag = 
A2~ block aibıcıdı T Ao_blockaabecods 
Aj0~block a1bıcıdı + A10—blockazbsceds 
ajo = A A (13) 
2—block a,b ic,d, T 2—blockagbe2ced2 
G= A2—block a1bıcıdı + Å2—blockazbzc2d2 
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Eq. (13) shows a complex set of high-order field harmon- 
ics and gradients. Based on previous research, it can be 
found that the high-order field harmonics a14 are small com- 
pared to ag and ajo. Therefore, the high-order field har- 
monics a14 are not considered for the time being, and the 
three physical quantities ag, a10, and gradient G are used as 
the objective function. Based on the characteristics of the 
model, this is a multiobjective optimization problem. Tra- 
ditional methods to deal with multiobjective optimization 
(MOO) problems usually convert such multiobjective prob- 
lems into single-objective optimizations by formulating some 
forms of cost functions or by converting some objectives into 
constraints. However, it can be extremely difficult to pro- 
vide a suitable weighting factor for determining preferences, 
and it is challenging to place proper constraints in engineer- 
ing practice. Several multiobjective evolutionary algorithms 
(MOEAs) have been proposed to address this issue [32]. 
In this study, the NSGA2 genetic algorithm (belonging to 
MOEAs) based on the pymoo platform in Python was used 
to solve the multiobjective optimization model. 

The calculation model of the genetic optimization algo- 
rithm is based on the biological evolution process of Darwin’s 
genetic selection and natural elimination and is a search algo- 
rithm that embodies biological genetics and the natural law 
of survival of the fittest. All individuals in the population are 
considered objects, and a coded parameter space is randomly 
and efficiently searched using selection, crossover, and muta- 
tion genetic operations. The initial population of the param- 
eter coding, fitness function, genetic operation, and control 
parameters constitute the core content of the genetic algo- 
rithm [33]. The focus of this study was the numerical calcu- 
lation of the magnetic field of the racetrack quadrupole coil. 
The genetic algorithm is not the focus of our research; it is 
only a means to obtain the optimal layout of the coil. There- 
fore, the Python database was used to directly call NSGA2, 
Problem, Minimize, and other function modules, and to mod- 
ify some parameters according to the problem to achieve our 
goal. 

In the solution model, the model characteristics were speci- 
fied before using the optimization algorithm to solve the prob- 
lem. The distance from the inside of the coil to the origin 
was 0.075 m, and the distance from the outside of the coil to 
the origin was 0.112 m. Two layers of racetrack coils were 
arranged between 0.075 and 0.112 m. The reference radius 
of the quadrupole coil was 0.050 m. These parameters are 
determined by the accelerator physics, magnet structure, and 
superconductors. As shown in Fig. 6, where each red ball 
represents a solution. It is not difficult to conclude that many 
noninferior solutions can be obtained using the genetic algo- 
rithm. 

In the two-dimensional image shown in Fig. 6, noninferior 
solutions were screened according to the requirements. There 
are two points at which the high-order field harmonics ag and 
aio Simultaneously satisfy the requirements of 5 x 1074 at 
the same time. Considering the influence of the gradient, the 
optimal solutions are listed in Table 1. 

Using the genetic algorithm, considering the limitations 
of space size, some solutions satisfying the conditions were 
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Fig. 6. Genetic algorithm solution results 


found. According to the solution results of the analytical cal- 
culation method, a quadrupole model with the same coil pa- 
rameters was established in the ROXIE software. A compari- 
son between the results of the analytical method and the Finite 
Element Method is also shown in Table 1. The results of the 
analytical and finite element methods are in good agreement. 
The values in parentheses in the FEM results are the simu- 
lation results of the corresponding normal quadrupole coils, 
which also prove the correctness of the relationship between 
the high-order field harmonics of the normal quadrupole coil 
and the skew quadrupole coil deduced above. 


By comparing the calculations of the analytical method 
with the finite element calculations, it was found that the cal- 
culation speed of the theoretical formula was faster. Further- 
more, when the finite element method searches for an optimal 
solution, the initial value, search range, and search step length 
of each parameter are required. This is the process of local 
optimization near the initial value. If the initial value is un- 
reasonable, the global optimal solution can easily be missed. 
However, the genetic algorithm is a process of global opti- 
mization, and it does not depend on the initial value [34, 35]; 
only the parameter range is required to find the optimal solu- 
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TABLE 1. Comparison of the analytical and FEM results. 


Parameter Analytical results FEM results 
a,=0.075 a2=0.094 
Block (m) b1=0.093 b2=0.112 


c1=0.0569394 
d1=0.0739649 


c2=0.02924927 
d2=0.07481863 


Current density (A/mm?) 1044 1044 
Gradient (T/m) -123.4 -123.4 
ag(10~*) -2.854 -2.867 (2.867) 
a1o(1074) -3.264 -3.263 (—3.263) 
a14 (1074) -0.118 -0.118 (0.118) 


tion. The genetic algorithm only requires a few seconds to de- 
termine the optimal solution. Therefore, the method of using 
the genetic algorithm to find the optimal solution is more reli- 
able and less time-consuming than the finite element method. 
Finally, according to the different application environments of 
the quadrupole coil, the genetic algorithm was used to obtain 
the optimal solution of the theoretical formula under different 
constraints. The theoretical formula-solving method is more 
flexible and convenient than the finite element method. 


HI. NUMERICAL CALCULATION ALGORITHM FOR 
THREE-DIMENSIONAL RACETRACK QUADRUPOLE 
COIL MAGNETIC FIELD 


The design of superconducting magnets can not only con- 
sider the cross-sectional magnetic field but also the magnetic 
field at the end of the magnet, which is very important and 
the most complicated and difficulty in the numerical calcula- 
tion of the three-dimensional magnetic field. The surface of 
the racetrack quadrupole coil is a regular plane and a cylinder 
that does not contain a complex spatial geometry. The end of 
the coil is illustrated in Fig. 7. The end coil can be regarded 
as the intersection of a plane and a cylinder. The equations 
for the plane and the cylinder are as follows: 


F (x,y,z) = x = Xo 


14 
G (x,y,z) =Y? + (z-L,)° = R? ea 


where F'(x,y,z) is the plane where the coil is located, 
G(x,y, z) is the cylindrical surface, Ls is half the length of 
the straight-line segment of the coil, and R is the bending ra- 
dius of the coil. 

The current direction in the conductor at the end of the coil 
is not parallel to the z-axis; therefore, to calculate the mag- 
netic field generated by the end coil, it is necessary to calcu- 
late the angle between the current and z-axis directions. 

The normal vectors of plane F and cylinder G are 


n= aa = (1,0,0) 

Ox’ Oy’ Oz (15) 
+ (0G 0G dG _ 
ng = (5 By? =) = (0, 2y, 2z — 2L,) 


Fig. 7. Intersection of three-dimensional racetrack coil end. 


30 The tangent vector of any point Q at the intersection AB is 


381 = ny X m = (0,2, — 2z, 2y) (16) 
The cosine of the angle between the current direction and 
the z-axis in the conductor can be obtained from the tangent 


vector. 
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As shown in Fig. 8, the coil ends are discretized. The x- 
axis is discretized into N parts, where each part is Ly = 
ba and the y-axis is discretized into P parts, where each 


N 3 
partis Lp = I£. M were dispersed along the z-axis. As the 
390 bending radius of the coil end is related to the number of parts 
1 dispersed along the y-axis, the length of the j-th part on the y- 
axis after dispersion along the z-axis is Lg = ce a 
To calculate the included angle between each end conduc- 
tor and the z-axis, the Q coordinates of the center point of the 
initial cross-section of the i-th conductor on the x-axis, the 
j-th conductor on the y-axis, and the k-th conductor on the 


z-axis are as follows: 
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The interface between 
the coil straight line 
segment and the coil end 


Each segment in 
the Z direction 


(d-c)/P 


Fig. 8. Discrete method of three-dimensional racetrack coil end. 


39 According to the multipole field component expressions of 
o each straight conductor, high-order field harmonics of the in- 
1 tegrated magnetic field generated by the i-th, j-th, and k-th 


conductor segments can be obtained [27]. 
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s4 For the straight coil section, the current in the conductor is 
5 parallel to the z-axis, that is, a=0; therefore, the above for- 
mula is also applicable to the straight section coil, and the to- 
7 tal integrated magnetic field of the racetrack coil is the super- 
ss position of the integrated magnetic fields generated by each 


s9 discrete conductor segment. 


4 


S 


406 


4 


So 


N P M 

Brn,coil = 5 5 5 (BL) (i jkjend + (BL) n straight 
i=1 j=1 k=1 
N P M 

An, coil = 5 5 5 (AL), (i,j,b)end F (AL)n straight 
i=1 j=1 k=1 


410 (20) 
Bn, coil and An coil are the n order normal and skew com- 
m2 ponents of the total integrated magnetic field, respectively. 
413 For normal and skew quadrupole coils, the relative magnetic 


44 fields b, and a,, can be described as follows: 
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According to the above analysis, the analytical expressions 
of the high-order field of the three-dimensional coils can be 
obtained. A three-dimensional coil layout with smaller high- 
order field harmonics of the integrated magnetic field can be 
determined by adjusting the length of the straight-line seg- 
ment. The length of the straight-line segment in the positive 
direction of the inner coil was 1075 mm, and the bending ra- 
dius of the end was the parameter c1; that is, the end shape 
was a regular semicircle. The length of the straight line in 
the positive direction of the outer coil was 1000 mm, and 
the bending radius of the end was parameter c2. The three- 
dimensional coil model is solved using the discrete sum- 
mation method. The number of discrete coils in the x-,y-, 
and z-directions affects the calculation accuracy. Combined 
with the performance of the computer and the accuracy of 
the model solution, the following discrete copies were deter- 
mined: The discrete numbers of inner and outer coils in the 
x-direction were 585 and 576, respectively, with a length of 
0.031 mm. The discrete numbers of the inner and outer coils 
in the y-direction were 553 and 1457, respectively, and the 
length of each part was 0.031 mm. The number of discrete 
parts of the inner coil along the z-axis was 370. Owing to the 
inconsistent lengths of the coil end, the minimum part length 
is 0.154 mm, the maximum part length is 0.2 mm, and the 
number of discrete parts of the outer coil on the z-axis is 187, 
the minimum part length is 0.156 mm, and the maximum part 
length is 0.4 mm. 


" 18.99T 


Fig. 9. Magnetic field analysis of OPERA three-dimensional coil 
model. 


The numerical results of the magnetic field calculation for 
the whole three-dimensional coil can be obtained by writing 
a code to calculate the discrete accumulation of the magnetic 
field. Racetrack quadrupole coils were constructed directly 
using the racetrack module in OPERA-3D and ROXIE. The 
high peak magnetic field at the end of the three-dimensional 
coil is caused by the multiturn cable winding, which can be 
reduced by grouping the cables and adjusting the length of 
the straight segments. Fig. 9 shows the calculation model of 
OPERA, and Table 2 compares the results of the analyti- 
cal algorithm, OPERA, and ROXIE for the high-order field 
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harmonics calculation. The difference between the 12-pole 
magnetic-field component results obtained using the analyt- 
ical formula and the software calculation results was three 
parts per million. The difference between the 20-pole mag- 
netic field component and 28-pole magnetic field components 
is small, which proves the feasibility of the three-dimensional 
magnetic field numerical calculation algorithm. 


TABLE 2. Calculation results of the three-dimensional algorithm 
for racetrack quadrupole coil. 


Field harmonics Analytical results OPERA ROXIE 


ag(10-*) 0.385 0.406 0.419 
a1o(1074) -3.060 -3.056 -3.056 
a14(1074) -0.100 -0.098 -0.099 


IV. PRELIMINARY ELECTROMAGNETIC DESIGN OF 


RACETRACK QUADRUPOLE MAGNETS 


Through a theoretical study of racetrack quadrupole coils, 
the geometric parameters of ideal racetrack coils satisfying 
physical requirements were obtained. The parameters a1, b1, 
C1, dı, a2, b2, C2, d2 for the skew racetrack quadrupole coil 
are listed in Table 1. Based on the conversion of the nor- 
mal and skew quadrupoles, the vertex coordinates of the two 
blocks in the skew quadrupole coil were rotated by 45° to be- 
come the vertex coordinates of the two blocks in the normal 
quadrupole coil. Rectangular superconducting cables are re- 
quired to complete the block filling. Because of the influence 
of the superconducting cable insulation layer, the coil param- 
eters must be fine-tuned during the electromagnetic design 
process. 

The electromagnetic design of the racetrack quadrupole 
magnet was performed using ROXIE software. Table 3 
lists the design requirements of a superconducting racetrack 
quadrupole magnet. 


TABLE 3. Design requirements of the large-aperture superconduct- 
ing quadrupole magnet. 


Parameters Value 
The inner radius of the quadrupole coils 75 mm 
Field gradient 132 T/m 
Reference radius 50 mm 
High order field harmonics <5x1074 


The double-layer normal racetrack coil structure is illus- 
trated in Fig. 10. The width of the cable was 18 mm, and 
the thickness was 1.5 mm. The inner layer of the coil had 
11 turns of superconducting cables and the outer layer had 28 
turns of superconducting cables. The peak field of the coil 
is 16.33 T. The calculated high-order field harmonics were 
within 3 x 1074. 

According to the particularity of the superconducting 
racetrack quadrupole coil, a bladder support structure was 
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Fig. 10. Calculation model of large-aperture and high-gradient race- 
track quadrupole coil 


adopted, as shown in Fig. 10. To adjust the magnetic 
field quality and form a reliable pre-stress structure, the 
cross-section of the iron yoke is divided, slotted, and perfo- 
rated [36, 37]. The simulation results are presented in Ta- 
ble 4. With the addition of an external iron yoke, the mag- 
netic field gradient increased to 132 T/m, and the peak field 
in the coil reached 17.23 T. The high-order field harmonics 
bg change in the positive direction. By adjusting the coil pa- 
rameters, the absolute value of the high-order field harmon- 
ics bg remained within 4 x 1074 after adding the iron yoke. 
When the magnetic field reaches 15-17 T, the current-carrying 
capacity of low-temperature superconductors is significantly 
reduced. Using a high-temperature superconductor is an ef- 
fective method of achieving high gradients [38, 39]. 


TABLE 4. Calculation results of large-aperture and high-gradient 
racetrack quadrupole coil. 


Parameters Coil Coil with iron 
Current density in strand 1166 A/mm? 

Field gradient 122.0 T/m 132.0 T/m 
Peak field in coil 16.33 T 17.23 T 
ag(10~*) 2.103 3.855 
a19(10~*) -2.825 -2.633 
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In the next example of a small-aperture racetrack magnet, a 
racetrack quadrupole coil was applied to complete the design 
of the superconducting quadrupole magnet QDO in the CEPC 
interaction region, and the CEPC adopted a double-ring col- 
lision mode [40]. There were two RF cavities [41, 42] and 
two interaction points on the main ring. Double-aperture con- 
ventional electromagnets were used in the main ring [43], and 
superconducting quadrupole magnets with high gradients and 
small apertures were used in the interaction region [44]. The 
design requirements of quadrupole magnet QDO are listed in 
Table 5. 


TABLE 5. Design requirements of small-aperture and high-gradient 
racetrack quadrupole coil. 


Parameters Value 
The inner radius of the quadrupole coils 20 mm 
Field gradient 136 T/m 
Reference radius 9.8 mm 
High order field harmonics <5x1074 


The cross-sectional layout of the coil obtained through pa- 
rameter optimization is shown in Fig. 11. The width of the 
cable was 3.5 mm and the thickness was 1.5 mm. Each pole 
had 12 turns of superconducting cable, and the current den- 
sity in the strands was 942.8 A /mm?. The high-order field 
harmonics are all within 3 x 10~*. To enhance the magnetic 
field gradient, we added an iron shell to the coil. The inner 
edge of the iron is a regular octagon. The radius of the inner 
tangent circle of the octagon is 35 mm, and the outer radius 
of the iron shell is 54 mm. 

Compared with the values of the quadrupole coils with- 
out an iron yoke in Table 6, when the iron yoke was added, 
the field harmonics increased slightly but were still within 
3 x 1074. The field gradient was 136 T/m. Supercon- 
ducting magnets in the interaction region are affected by 
radiation from the particle beam. Considering the stabil- 
ity of the superconducting magnet, a small-aperture, high- 
gradient quadrupole magnet can also be designed with a high- 
temperature superconductor to leave an ample margin. 


TABLE 6. Calculation results of large-aperture and high-gradient 
racetrack quadrupole coil. 


Parameters coil coil with iron 
Current density in strand 942.8 A/mm?” 

Field gradient 102 T/m 136 T/m 
Peak field in coil 3.328 T 4.073 T 
ag(10~*) -0.572 -2.018 
a19(10~*) -1.762 -1.312 


These two racetrack quadrupole magnet examples are the- 
oretically feasible. In our study, only two rectangular blocks 
were used to adjust the high-order field harmonics, and the 
systematic field harmonics achieved a precision of 1074. In 
specific engineering applications, more rectangular blocks 
may be required to satisfy the design requirements. In sum- 
mary, using a superconducting cable to fill the block area 
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Fig. 11. Calculation model of small-aperture and high-gradient race- 
track quadrupole coil 


54 in the large-aperture high-gradient model and small-aperture 
sat high-gradient model can satisfy our design requirements, and 
sa2 the magnetic field quality accuracy reaches the order of 10~4. 
sss According to the superconductor critical current density 


a 
= 


a 
= 


580 
581 
582 


583 


10 


curve, the peak magnetic field and current-carrying density 
in the racetrack coil with a large aperture and high gradient 
exceed the critical performance of low-temperature supercon- 
ductors NbTi and Nb3Sn, and the isotropic high-temperature 
superconductor Bi-2212 is the best choice [45, 46]. Racetrack 
coils with small apertures and high gradients can simultane- 
ously use both low- and high-temperature superconductors si- 
multaneously. The advantages of using a low-temperature su- 
perconductor are that the coil structure is simple and conve- 
nient to manufacture and the operating margin of the magnet 
is relatively sufficient. The use of a high-temperature super- 
conductor can further increase the operating margin, stability, 
and anti-interference ability of magnets. 
V. SUMMARY 


Through research on a superconducting racetrack 
quadrupole coil, analytical formulas for high-order field 
harmonics and field gradients with respect to the current den- 
sity and coil, geometric parameters were obtained, and the 
analytical numerical calculation results of the magnetic field 
were consistent with the finite element calculation results. 
A genetic algorithm was applied to obtain a solution for the 
coil geometry parameters with field harmonics on the order 
of 1074. Finally, preliminary electromagnetic designs of 
large-aperture and small-aperture racetrack superconducting 
quadrupole magnets are presented. 

The racetrack superconducting quadrupole coil had a sim- 
ple shape and manufacturing process, in which both high- 
and low-temperature superconductors can be used. A race- 
track quadrupole coil with a simple geometric structure is 
beneficial for strain-sensitive high-temperature superconduc- 
tors and is an effective way to realize ultra-high magnetic field 
gradients. For low-field magnets, the use of low-temperature 
superconductor racetrack coils can simplify the process and 
reduce costs. High-temperature superconductors can also be 
used in low-field racetrack magnets to increase the magnet 
operating margin and improve magnet stability. Racetrack 
superconducting quadrupole magnets are a frontier research 
direction, and the research presented in this paper provides 
a new concept for the design of superconducting quadrupole 
magnets. 


584 VI. APPENDIX 


ses Each order multipole field component A,, of the entire racetrack quadrupole coil was obtained. The analytical expressions 
ses for the most important multipole field components, Ag and Ajo are 
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